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ABSTRACT 

We propose to determine the stellar velocity dispersions of globular clusters in the 
outer halo of the Milky Way in order to decide whether the dynamics of the universe 
on large scales is governed by dark matter or modified Newtonian dynamics (MOND). 
We show that for a number of Galactic globular clusters, both the internal and the 
external accelerations are significantly below the critical acceleration parameter ao 
of MOND. This leads to velocity dispersions in case of MOND which exceed their 
Newtonian counterparts by up to a factor of 3, providing a stringent test for MOND. 
Alternatively, in case high velocity dispersions are found, these would provide the first 
evidence that globular clusters are dark matter dominated. 
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1 INTRODUCTION 

Dark matter is now generally believed to be the domi- 
nating mass compon e nt of the universe, starting with the 
discovery of IZwickvl dl93St) that the speed of galaxies in 
the Coma cluster is too large to keep them gravitation- 
ally bound unless they are much heavier than one would 
estimate on the basis of visible matter alone. Although the 
currently favoured ACDM mode l has proven to be r emark- 
ably successful on large scales dSpergel et alJl2005t . high- 
resolution N-body simulations are still in contradiction with 
observations on subgalactic scales where they predict or- 
ders of magnitude more substructure than what is seen 
jMoore et alJ ll99St iKlypin et all Il999l) and also a wrong 
spati al distribution of the subhalos llKroupa. Theis fc Boilvl 
l2004h . Additional arguments regarding supporting and con- 
tradictory observation s of structure formation in CDM mod - 
els are presented in iGreb el. Gallagher fc Harbeckl (120051) 
and IGrebel fc Gallagherl j2004l) . The above discrepancies 
might be resolved if more realistic simulations that can bet- 
ter treat the dynamics of the interstellar gas and feedback 
processes become available, or they could show that our cur- 
rent understanding of cosmology and large scale structure 
formation is still missing important ingredients. 

An alternative to the dark matter hypothesis could be 
Modified Newtonian Dynam ics (M OND), which was pro- 
posed by iMilgroml lll983 a b) and Bekens tein fc Milgroml 



( 1984) as a way to explain the rotation curves of galaxies 
without the need to assume large amounts of otherwise un- 
seen dark mass in the outer parts of galaxies. According to 
MOND, Newtonian dynamics breaks down in the limit of 
very weak accelerations, and the acceleration a experienced 
by a particle is given by the following (heuristic) equation: 



(1) 



M — a = 9n 

a 



where 



H(x) 



x if x « 1 
1 if x » 1 



(2) 



* e-mail: holger@astro.uni-bonn.de (HB); grebel@astro.unibas.ch 
(EKG); pavel@astro.uni-bonn.de (PK) 



Here gh is the standard Newtonian acceleration and ao a 
constant , observationally determined to be ap = 1.2 ■ 10~ 8 
cm/sec 2 jBegeman et al]ll99ll ISanders fc McQaugbll2002f ). 
The above formulation of MOND conserves angular momen- 
tum and energy only for spherical mass configurations, but 
a more general formulati on exists which obeys these c onser- 
vation laws for all cases feekenstein fc Milgromlfl984h . 

Predictions from MOND have been shown to be in good 
agreement with the observed rotati o n curves of galaxies 
jBegeman et alJ Il99ll : ISandersl I199& ISanders fc Verheiienl 
1998) and can also explain the velocities of galaxies in 
groups wi th reasonable M/L values for individual galaxies 
jMilgromlf2002l) . In addition, IMilgroml <1995ft showed that 
the velocity dispersion of 7 dwarf galaxies which was avail- 
able at the time was compatible wi th the predic tions from 
MOND, a conclusion also found bv iLokasl (2002) with up- 
dated data for the Fornax and Draco dwarf galaxies. These 
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successes are remarkable given the fact that unlike the dark 
matter hypothesis, MOND has only one free parameter (ao) 
that can be adju sted to explain observations. In addition, 
iBekensteinl <l2004f) has recently developed a relativistic for- 
mulation of MOND, putting the theory on a more solid the- 
oretical basis. 

Part of the trouble in deciding whether MOND or dark 
matter is the better candidate in explaining the velocities of 
stars on galactic scales stems from the fact that cosmolog- 
ical structure formation is still not sufficiently understood, 
so within certain limits, the mass and size distribution of 
dark matter can be adjusted to fit the observational data. 
It is therefore highly desirable to test MOND for objects in 
which no dark matter is thought to exist. An ideal candi- 
date for such objects are globular clusters, which still form 
today as a result of collis ions between gas clouds durin g ma- 
jor mergers of galaxies iWhitmore fc Schweizerl|l995f) . and 
which are beli eved to have formed in the same way in the 
early universe (lAshman fc Zepflll992D . 

Testing gravity for low accelarati ons by using nearb y 
globular clusters was already tried bv lScarpa et alJ {2003). 
However, nearby clusters experience an acceleration from 
the Milky Way that is larger than the critical MOND con- 
stant ao and should therefore by governed by Newtownian 
dynamics if MOND is correct. 

In the present paper we investigate the effects of MOND 
for a number of low-mass globular clusters in the outer halo 
of the Milky Way. It will be shown that if MOND is true, 
these clusters would have mass-to- light ratios far larger than 
their Newtonian values, allowing an independent test of the 
predictions of MOND. The paper is organised as follows: 
In section 2 we calculate expected velocity dispersions and 
mass-to-light ratios for a number of globular clusters for 
the Newtonian case and in the MOND regime. Section 3 
discusses how the validity of MOND could be constrained 
if the velocity dispersions of these clusters would be known 
and in section 4 we will draw our conclusions. 



2 EVALUATING THE EFFECT OF MOND 
FOR GLOBULAR CLUSTERS 

2.1 LOS velocity dispersions for different cases 

In case of Newtonian gravity, the virial theorem connects 
the mass Mc, radius r and average velo city dispersion a of a 
clust er through the following equation iBinnev fc Tremaind 
Il986ft . eq. 4-80a: 



GM C 



(3) 



where r v is the virial radius, which in many stellar systems 
can be approximated by the three-dimensional half-mass ra- 
dius Th as r v « 2.5rj, if the clusters are stationary systems 
and sufficiently unperturbed that the assumption of virial 
equilibrium is valid. A similar relation exists between the 
three-dimensional half-mass radius and the easier to observe 
two-dimensional, projected half-mass radius rh v - Th v = 7TV1 
with 7 » 0.74. 

If we assume an isotropic velocity distribution, the 
line-of-sight velocity dispersion is related to the three- 
dimensional one through ulos = cr/y3, and we obtain the 
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Figure 1. Ratio of the line-of-sight velocity dispersion of stars in 
different King models with the prediction from eq. 21 (solid line). 
For all models, the velocity dispersion agrees with the prediction 
to within 5%. 



following equation for the observed velocity dispersion of a 
star cluster in case of Newtonian gravity 



&los,n = 0.335 



GM C 



(4) 



In order to test t he accuracy of eq. 0] we created N- 
body representations of lKind (1966) models with dimension- 
less central concentrations in the range 3 ^ Wo ^ 15 and 
checked our theoretical estimate against the JV-body data. 
King models with concentrations in this range are usually 
used to represent density profiles of globular clusters. Fig. 1 
compares the average velocity dispersion of stars in these 
models with the prediction of eq. 01 For all models, the 
agreement is within 5%, which is accurate enough for our 
purpose of predicting velocity dispersions in real clusters. 

In case of MOND, solutions exist only for several special 
cases since the function fi(x) is defined only for limiting val- 
ues of x. If the acceleration a that stars experience is much 
larger than ao, fj, « 1 and the MOND solution is the same as 
in the Newtonian case. For this to be true it does not mat- 
ter if the acceleration a is the internal acceleration m n t due 
to stars in a cluster or if it comes from an external grav- 
itational field a e xt- External accelerations are for example 
important for all globular clusters which have galactocen- 
tric distances Rgc < few • 10 kpc or for the dynamics of the 
solar system. In both cases the external acceleration due to 
the Milky Way is larger than ao and the dynamics can be 
described by Newtonian gravity, no matter how small the 
internal accelerations are. Hence, MOND can be important 
only for star clusters in the outer halo of the Milky Way. 

If both a,i n t and much smaller than ao one 

is in the deep MOND regime with fi(x) w x. If in this 
case the external acceleration is larger than the internal one 
a,ext » ai n t, the system is again nearly Newtonian but with 
an effective gravitational constant G that is larger than the 
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standard Newtonian one by a factor ao / a ex t iMi lgromll986l) . 
The line-of-sight velocity dispersion is therefore equal to 

nw \ gm c fW , . 

°los,mi — 0.665* . . (5) 

V r hp y a ext 

Finally, if a ext « a» n t and both are smaller than ao, 
the cluster is effectively isolated and the acceleration of the 
cluster stars is given by gu = ^/aogjv where is the ac- 
celeration in the Newtonian case. The line- of-sight velocit y 
dispersion of a star cluster is then given bv dMilgroml lioQl: 

olosmi = 0.471 (a GM c ) 1/4 , (6) 
i.e. independent of the clusters radius. 



a ext = (7) 

and hold in both the Newtonian and the MOND case. The 
line-of-sight velocity dispersions for the MOND case were 
calculated from eqs. 5 and 6, depending on whether the in- 
ternal or external acceleration is larger. While this proce- 
dure is not strictly valid within the framework of MOND, it 
is accurate enough to give an estimate of the effect MOND 
would have. The M/L values in the last column finally are 
calculated from the olos,m values assuming that an ob- 
server interpretes them for the Newtonian case. They should 
be compared with the Newtonian input value of M/L — 2 
which we assumed. 



2.2 Application to globular clusters 

Table 1 shows the predicted velocity dispersions of globular 
clusters for the Newtonian and the MOND case. We have 
calculate d the dispersio ns for all galactic globular clusters in 
the list of lHarrisI (119961) . but present only those in Table 1 for 
which there is a noticeable difference between the two cases. 
Clusters in Table 1 are generally far away in the galactic halo 
so that the external acceleration due to the Milky Way is 
small, and also have small masses and large half-mass radii 
so that their internal accelerations become small. 

Galactocentric distance s Rac and projected half-mass 
radii are taken from lHarrisI (|l99(f). We have assumed that 
for the clusters of Table 1 the half-light radius is equal to 
the half-mass radius i.e. mass follows light. This is probably 
a good assumption since the half-mass relaxation times of 
most clusters in Table 1 are of the order of a Hubble time or 
even larger, so dynamical evolution is not likely to play an 
important role for these clusters. Cluster masses were cal- 
culated from the absolute luminosities by assuming a stel- 
lar mass-to-light ratio of M/L = 2. This value is in agree- 
ment with observed mass-to-light r atios o f globular clusters 
llMandu shev. Spassova fe Staneva l Il99lt iPryor fc Mevlanl 
ll993h . lM^ndushev^p^iaso'v^^'taneval Jl99lf) for example 
fitted single-mass King models to a sample of 32 clusters 
with reliable central velocity dispersions a nd determined an 
averag e M/L — 1.21 for their sample, while IPrvor fc Mevlanl 
(1993) used multi-component King-Michie models for 56 
clusters and estimated an average global mass-to-light ra- 
tio of M/L — 2.3. Since multi-component models can take 
the effect of mass segregation into account, the latter value 
is probably closer to the truth. 

The tidal radii Rt in Table 1 were calculated from 
the cluster masses by assuming a flat rotation curve of 
the Milky Way with Va = 200 km/sec for all galacto- 
centric distances. For most clusters they are a factor of 5 
to 10 larger than the half-mass radii, which means that 
tidal effects play no significant role for the internal dy- 
namics of the clusters in our sample. This is a signif- 
icant improvement compared to the situation for dwarf 
galaxies which generally have larger Rh/Rt ratios and 
for which there is an ongoing discussion to which ex- 
tent the high-M/L ratios found are caused by tidal effects 
(Kroupa 1997; Klcsscn fc Krouj3a|^)9S|; |Odgriki rchen et alJ 
l200ll:lFleck fc Kuhj2003UWilkinson et al.ll2004) ~ 

Internal accelerations were calculated at the clusters 
half-mass radii from gM = y'aogjv, assuming that MOND is 
true, while the external accelerations were calculated from 



3 DISCUSSION 

From Table 1, it is evident that if MOND is true, olos 
and the deduced M/L are increased by a significant amount 
over the Newtonian case. Three clusters would have M/L 
values in excess of 10, which is far larger than what is 
fo und in any galactic globular cluster. In the analysis 
of iMandushev. Spassova fc Staneval dl99lf) for example no 
cluster had an M/L larger than 3.0 and low-mass clusters 
with IgMc < 4.5 generally h a d M/ L < 1.5. Similarly, in the 
analysis of IPrvor fc Mevlanl (1993) low-mass clusters with 
IgMc < 5 had an average M/L = 1.8 and none had an 
M/L larger than 3.2. If a high M/L ratio would be found 
in any cluster, it could therefore be interpreted in only one 
of three ways: 

(i) The initial mass-function of this cluster was radically 
different from ordinary GCs and heavily weighted either to- 
wards high-mass stars which have by now turned into com- 
pact remnants, or towards very low-mass stars in order to 
create such a large M/L. 

(ii) The cluster contains a significant amount of cold dark 
matter. 

(iii) Newtonian dynamics has to be modified in the limit 
of low accelerations. 

Possibility (i) can almost certainly be excluded since such 
high M / L ratios have not be en found for any s tellar 
population to da te. In addition, Ste tson et alJ lll999f) and 
ISaraiediiil lll997l) determined ages and metallicities for four 
of the clusters in our list (Eridanus, Pal 3, 4 and 14) which 
place them among other galactic globular clusters, indicat- 
ing that the r e is n othing unusual about these clusters. 

iPeeblesI dl984f) and others after him have suggested that 
globular clusters formed in cold dark matter mini-halos in 
the early universe, in which case they should contain sig- 
nificant amounts of dark matter. So possibility (ii) seems 
entirely possible and a detection of dark matter dominated 
globular clusters would be a direct confirmation of this idea. 
In this case halo globular clusters could help to bridge the 
gap between the predicted number of DM sub-clumps in the 
halos of galaxies and th e observed numb er of clumps in the 
form of dwarf galaxies (ICote et al.ll2002T) . Since the forma- 
tion history and dark matter content would basically be the 
same, halo globular clusters could in such a case be viewed 
as smaller-sized versions of dwarf galaxies. 

However, observations of interacting and starburst 
galaxies have shown that the formation of globular clus- 
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Table 1. Distant globular clusters for which the predictions of MOND and Newtonian dynamics differ. For explanation of symbols see 
text. 



Cluster 


lg M c 


rh P 


Rrp 


Rac 




&ext 


&LOS,N 


°LOS,M 


M/L 


Name 


[Mq] 


[PC] 


[PC] 


[kpc] 


[10 -9 cm/s 2 ] 


[10 -9 cm/s 2 ] 


[km / sec] 


[km/sec] 


MOND 


AM 1 


4.10 


17.7 


151.3 


123.2 


1.84 


1.05 


0.58 


1.77 


18.2 


Eridanus 


4.27 


10.5 


145.5 


95.2 


3.79 


1.36 


0.93 


1.96 


8.8 


Pyxis 


4.52 


15.6 


101.1 


41.7 


3.38 


3.11 


1.01 


2.25 


9.9 


Pal 3 


4.50 


17.8 


173.5 


95.9 


2.90 


1.35 


0.92 


2.23 


11.6 


Pal 4 


4.63 


17.2 


212.0 


111.8 


3.48 


1.16 


1.09 


2.40 


9.6 


Pal 14 


4.11 


24.7 


103.4 


69.0 


1.33 


1.88 


0.50 


1.27 


12.8 


Pal 15 


4.42 


15.7 


87.6 


37.9 


2.98 


3.42 


0.90 


1.68 


7.0 


Arp 2 


4.34 


15.9 


56.3 


21.4 


2.68 


6.06 


0.81 


1.14 


4.0 



ter is still an ong oing process, triggered main ly by ma- 
jor galaxy mergers iWhitmore fc Schweizerll995t) . Similarly, 
the relatively high m etal licities and young er ages found by 
IStetson et alJ lll999li and ISaraiedinl Jl997t) speak against a 
primordial formation sce nario of these globular clu sters and 
for an accretion scenario flMackev fc Qilmorel2004h . It there- 
fore appears unlikely that the clusters in Table 1 contain a 
significant amount of dark matter. 

Another possible interpretation of high velocities would 
be that Newtons law of gravity has to be changed for 
low accelerations. Observations of additional halo clusters 
for which MOND does not predict a velocity enhancement 
would help to distinguish between the dark matter hypoth- 
esis and MOND. Also, checking the evolutionary state of 
the clusters which are in the MOND regime could help de- 
cide between MOND and dark matter since the relaxation 
time in the MOND regime is significantly smaller than in 
the Newtonian case ICiotti fc Binnevll2004h . meaning that 
the clusters would be dynamically more evolved if MOND 
is true. 

A second possible outcome is that the observed veloc- 
ity dispersion is in complete agreement with the Newtonian 
value and far lower than what MOND predicts. In such a 
case MOND would probably have to be discarded in its 
present form since it would be impossible to reconcile the 
high rotational velocities and velocity dispersions seen in 
galaxies with the low velocities seen in the clu s ters o f Ta- 
ble 1: Assuming Newt onian gravity iKochanekl (119961) and 
ISakamoto et al. ll2003f) found from an analysis of the or- 
bital motion of galactic satellites that the total mass of the 
Milky Way within 50 kpc is respectively (4.9 ± 1.1) ■ IO^Mq 
and (5.5 ±0.2) • lO n M , giving rise to a rotational velocity 
of Vg = 210 km/sec. On the other hand, the mass of the 
Milky Way in the disc and bulge in the form of visib le mat- 
ter is of order 5.5 • 10 10 M Q JPehnen fc Binne\ll998h . giving 

/ \ 1/4 

rise to a rotational velocity of Vg = 164 I 10 -&cm/ sec' 1 ) 

km/sec within the framework of MOND (eq. |SJ. Values of 
ao below 5 • 10 -9 cm/sec 2 seem therefore to be very hard 
to reconcile with the orbital velocities of galactic satellites, 
even given the uncertain ties in the above numb ers. Similarly 
the galaxies analysed bv lBegeman et al] (1199 ll) give a mean 
ao = (1.21 ± 0.27) ■ 10 -8 cm/sec 2 , excluding values of ao 
below 6.7- 10 -9 cm/sec 2 at the 2a level. On the other hand, 
even an ao as low as 5- 10 -9 cm/sec 2 would still give a veloc- 
ity dispersion of 1.42 km/sec for AM 1 and 1.80 km/sec for 
Pal 3, which is twice the Newtonian value. Hence, velocity 



dispersion measurements will have the power to falsify or 
support MOND. 



4 CONCLUSIONS 

We have calculated the velocity dispersions for a number of 
globular clusters in the halo of the Milky Way and shown 
that in case of MOND they would significantly exceed the 
corresponding Newtonian values, allowing a test of MOND 
and dark matter theories for a new class of objects and for 
length scales one order of magnitude smaller than where 
they could be tested before. We showed that interesting in- 
sights can be obtained about the formation of globular clus- 
ters and the role of dark matter almost independently of the 
actual results, so an observational effort to determine the ve- 
locity dispersions for clusters in our sample would be highly 
rewarding. Measurements of this kind are feasible with the 
current generation of 8 to 10-m class telescopes. 

Compared to dwarf galaxies, the studied globular clus- 
ters also have half-mass radii one order of magnitude smaller 
than their tidal radii, which means that they are effectively 
isolated. Their velocity dispersions are therefore much more 
straightforward to interpret since tidal effects are not likely 
to play an important role. Also, unlike tests for MOND 
based on the velocity dispersion of stars in the halos of 
nearby globular clusters, the interaction of the cluster stars 
with the galactic tidal field and the proper separation of 
members and non-members is not likely to be a problem for 
our clusters. 
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